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1 .  INTRODUCTION 

Suppose  we  have  a  connected  network  G(N,A)  consisting  of  a 
set  of  nodes  N  and  a  set  of  arcs  A.  Let  the  integers  i  -  1,2,..., P 
represent  the  nodes  and  the  two-tuples  (i,j)  (i  -  1,2, . . .  ,P  ; 
j  «  1,2,...,P;  and  i  i  j)  represent  the  ercs.  Let  node  1  corre¬ 
spond  to  the  source  and  node  n  correspond  to  the  sink.  The  arcs 
are  assumed  directed  so  that  the  order  (l,j)  implies  an  arc  directed 
from  node  i  to  node  j . 

Let  ^  0  represent  the  initial  flow  capacity  of  arc  (l,j), 
0  represent  the  added  flow  capacity,  and  X^  i  0  represent 
the  actual  flow  in  the  arc.  Let  Q  represent  the  net  flow  through 

P 

the  network  from  node  1  to  node  Let  B  i  0  be  the  total 

resource  budget  available  and  let  a^  i  0  be  the  cost  per  additional 
unit  of  capacity  added  to  (i,J). 

The  capacity  expansion  problem  for  a  flow  network  was  first 
studied  by  Fulkerson  (1959)  .  He  referred  to  it  as  the  parametric 
budget  problem  and  stated  it  as  follows:  Find  nonnegative  values 
of  and  which 


maximize  Q, 

(la) 

subject  to  £ 

J 

(W " Q’ 

l 

j 

(Xij_Xji)  ■  0  (i  -  2,3 . P-D, 

(lb) 

I 

l 

J 

(Xpj-Xjp)  ■  “  Q. 

1 

X  -Y  *M  ((i,  j)  €  A) , 

(lc) 

l 

^A  aijYij  *  B> 

(Id) 

/ 


a 


2 


In  contrast  to  Fulkerson  we  write  the  budget  constraint  (Id)  in 
inequality  form  although  it  is  obvious  that  equality  will  hold  at 
optimality  if  Y  is  not  required  to  be  integer  valued. 

An  easy  way  to  solve  (1)  is  to  construct  a  network  having 
two  parallel  arcs  in  place  of  each  original;  one  with  a  cost  of 


zero  and  one  with  a  cost  of  aij  ‘  The  capacity  of  the  arc  having 
zero  cost  would  be  while  the  other  would  be  infinite.  A  minor 

modification  of  tne  Primal-Dual  algorithm  of  Ford  and  Fulkerson  (1957) 


can  then  be  used  to  find  the  maximal  flow  which  can  be  allocated  for 


the  available  budget.  The  resulting  algorithm  is,  in  fact,  the  one 
given  by  Fulkerson  (1959)  in  his  presentation  of  the  parametric  budget 
problem. 

In  this  paper  we  would  like  to  first  consider  the  special  form 
of  (1)  when  -  0.  We  will  state  the  associated  optimal  solution 

and  prove  it  constructively  using  the  decomposition  principle  for 
linear  programs  (Dantzig  (1963)).  We  will  then  present  an  algorithm 
for  bolving  the  general  form  of  (1)  which  makes  use  of  the  topological 


dual  of  the  network. 


3 


2.  OPTIMAL  SOLUTIONS  WHEN  M  -  0 

The  form  of  (lc)  reduces  to 

XU  *  YU 

when  -  0  and,  in  fact,  is  equality  for  all  arcs  since  the 

budget  will  be  used  up  in  obtaining  the  maximal  possible  flow.  We 
can  'herefore,  drop  (lc)  from  further  consideration  if  we  rewrite 
(Id)  as 

I  IA  .tjXy  .  B.  (2) 

The  optimal  solution  to  (la,b)  and  (2)  is  given  by  the 
following  theorem. 

Theorem:  The  maximal  klow  will  be  bent  eve a  that  chain  C 

a 

dixected  hnom  &ou/ice.  to  sink  which  is  the  bhoxtest  fjoute  when 
the  axe  length  o<$  (i,j)  is  bet  equal  to  a^.  The  value  o 6 

the  maximal  fa low  and  the  ^’.ou>  capacity  in  each  axe  o&  that 
chain  is  given  by 

B 

max  Q  -  -  -  Y  . 

I  k  •*,  11 

m 

The  &law  capacities  in  all  othex  axes  axe  zexo. 

Proof:  Let  us  begin  by  adding  an  arc  (P,l)  to  the  network 
and  then  write  (lb)  as  EX  -  0  where  E  is  the  node-arc  Incidence 
matrix  of  our  augmented  network,  X  is  the  augmented  column  vector 
of  arc  flows  (with  Xp^  in  the  last  position),  and  0  is  null 


A 


column  vector.  We  write  (2)  as  aX  £  B  where  a  ■  (a^)  a  row 


vector  with  ap^  ■  0  in  the  last  position. 

find  X  ^  0  which 

Our  problem  is  then  to 

maximize  Xp^, 

(3a) 

subject  to  EX  ■  0, 

(3b) 

and  aX  £  B. 

(3c) 

In  decomposed  form  we  will  combine  (3a)  and  (3c)  into  the 

m.  ter  program.  The  flow  conservation  equations  (3b)  and  X  i  0 

will  be  the  constraint  set  of  the  subprogram.  The  master  program 

is  therefore  to  fin.’  A.  i.  0,  x  i  0  which 

i  s 


maximize  V  bX^A,, 
i  1 

(Aa) 

subject  to  £  aX^A^  +  xg  ■  B. 

(Ab) 

If  the  network  has  m  arcs  then  b  is  an  1  *  m  +  1  vector  con¬ 
sisting  of  zeros  in  the  first  m  positions  and  a  one  in  the  last 
position.  The  subprogram  is  to  find  X  i  0  which 

maximizes  [b  -  ira]X,  (5a) 

subject  to  EX  ■  0. 

We  are  using  tt  to  represent  the  simplex  multiplier  associated  with 
(Ab) .  Each  X^  us  id  in  the  master  program  is  a  homogeneous  solu- 
to  the  subprogram  and  therefore  we  have  no  £  A.  ■  1  constraint  in 


the  master. 
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The  iterative  solution  procedure  begins  by  finding  the  basic 

feasible  solution  x  ■  B  to  the  first  restricted  master  which  has 

s 

only  one  variable;  namely,  xg.  The  multiplier  tt  is  then  determined 
and  used  to  provide  the  first  objective  function  for  the  subprogram. 

If  a  solution  to  the  subprogram  can  be  found  such  that 

[b  -  wa]X  >  0  (6) 

then  we  use  that  solution,  call  it  X^  ,  to  generate  a  new  coeffi¬ 
cient  vector  for  the  restricted  master.  The  second  restricted  master 
program  would  then  be:  Find  A^  i  0,  xg  2  0  which 

maximize  bX^A^, 

subject  to  aX^A^  +  xg  ■  B. 

A  new  value  of  n  would  be  determiced  and  used  to  generate  a  new 

objective  function  for  the  subprogram.  A  new  solution  to  the  subproblem, 

(2) 

X  ,  would  be  determined  and  tested  using  (6) .  If  it  fails  the 

test  then  the  solution  process  terminates.  If  it  passes  the  test 

then  it  Is  used  to  create  a  new  restricted  master  program  involving 

variables  A,,  A„,  and  x  . 

1*  2*  s 

The  process  terminates  after  the  Jc^1  iteration  where  X^ 
is  the  first  solution  from  the  subprogram  to  fail  the  test.  The 
optimal  solution  X*  can  be  obtained  from  (7) . 

X*  -  l  AX(i)' 
i-1 


(7) 


6 


The  basic  feasible  solution  to  the  first  restricted  master 

was  observed  to  be  x  -  B.  Because  there  is  a  zero  coefficient 

s 

for  xg  in  (4a)  we  get  it  *  0  as  the  value  of  the  simplex  multi¬ 
plier  for  the  restricted  master  and  the  objective  function  for  the 
subproblem  is  therefore  bX  ■>  X^.  A  homogeneous  solution  to  the 
first  subprogram  is  easily  provided  by  selecting  any  cycle  directed 
from  the  source  to  the  sink  and  back  to  the  source  which  includes 
arc  (P,l)  and  allows  as  much  flow  as  possible  to  pass  over  it. 
Since  there  are  no  flow  capacity  restrictions  in  the  subproblem  we 
could  have  X^  -*•  °°.  If  the  solution  to  (3a,  b,  c)  is  unbounded 
this  condition  will  be  provided  by  X^  -*•  °°  for  some  solution  in 
(4a,  b)  so  we  need  only  send  1  unit  over  the  cycle  in  the  sub¬ 
problem  to  indicate  the  homogeneous  solution  route.  Thus,  we  return 
to  the  master  with  the  vector  X^  consisting  of  plus  ones  corre¬ 
sponding  to  arcs  on  the  cycle  selected  and  zeros  elsewhere.  Note 
that  [b  -  ira]X  =  1  >  0  so  (6)  is  satisfied. 

The  next  restricted  master  program  is:  Find  X..  0,  x  2:  U 

which 

maximizes  X^ , 

subject  to  aX^X^  +  xg  ■  B. 

The  optimal  basic  feasible  solution  to  this  problem  is  obviously 

x  =  0  and 
s 


7 


where  denotes  the  cycle  just  chosen  in  the  subprograr.  The 

associated  value  of  it  is  given  by  (8)  . 


I  l 


(8) 


ij 


The  objective  of  the  subprogram  is  then  to  maximize 


X 


PI 


aX 


(9) 


In  other  woras,  the  "cost"  for  sending  flow  over  an  arc  in  the  sub¬ 
program  is 


l  l 


C 


1 


(10) 


for  all  original  arcs  of  the  network  and  -1  for  arc  (P,l). 

The  maximum  value  of  (9)  for  a  flow  of  one  unit  would  be 
obtained  by  selecting  that  cycle  which  had  the  cheapest  "cost"  from 
source  to  sink  and  which  returns  to  the  source  via  arc  (P,l). 
Suppose  we  denote  this  cycle  by  C The  associated  solution  will 
be  X^  if  (6)  is  satisfied.  And  (6)  will  be  satisfied  if  (9)  is 
positive.  The  value  of  (9)  for  this  solution  is 


1  - 


l  l 


‘ij 


l  l 


ij 


» 


8 


and  therefore  (9)  will  be  positive  only  if 


(ID 


Thus  we  introduce  a  new  coefficient  vector  into  the  restricted 
master  only  if  the  route  we  found  on  the  second  visit  to  the  sub¬ 
program  was  shorter  in  the  sense  of  the  a^'s  than  on  the  first 
visit. 

Suppose  we  assume  that  (11)  is  true.  The  new  restricted 
master  then  has  the  form:  Find  A.  i  0,  A-  £  0,  x  2  0  which 

X  Z  8 

maximizes  A^  + 

subject  to  aX^A^  +  aX^A£  +  xfl  “  B. 

The  optimal  solution  to  this  master  is  A^  *  xg  ■  0  and 


n 


becavse  of  (11).  The  associated  value  of  tt  is 


7t 


9 


and  the  new  subprogram  objective  function  is 


*P1- 


aX 


9 


The  new  "cost"  of  sending  a  flow  over  an  arc  is  therefore 


l  l 


C 


2 


(12) 


for  arcs  of  the  original  network. 

When  we  compare  (12)  with  (10)  we  see  that  only  the  denomina¬ 
tors  differ.  Recalling  that  a^  i  0  for  all  arcs  we  then  realize 

(2) 

that  the  next  solution  to  the  subprogram  will  be  identical  to  X 

because  we  will  get  the  same  cheapest  route  cycle.  The  inequality 

(6)  will  be  violated  and  the  solution  process  will  terminate.  The 

cycle  C_  therefore  consists  of  the  chain  C  defined  in  the 

2.  m 

theorem  and  arc  (P,l).  The  value  of  ^  is  therefore 


x - 2 -  (13) 

2  II  « 

C 

m 


because  ap^  ■  0. 

(2) 

The  optimal  solution  to  (3a,  b,  c)  is  X*  -  where  ^ 

is  given  by  (13) .  Thus  the  arcs  on  the  chain  having  the  shortest 
length  in  the  sense  of  a^'s  between  the  source  and  sink  should 
have  a  flow  of  and  their  increased  arc  capacities  should  be  X^. 

The  maximal  flow  under  these  conditions  will  be  X ^  through  the 
network. 

The  following  corollary  is  an  immediate  consequence  of  the 
theorem  if  a^  ■  k  for  all  arcs. 


10 


Conollaj iy:  l  $  =  k  60fl  all  ones  06  the.  network  then  the 

optimal  solution  Is  to  spend  all  o£  the  budget  on  the  chain 
having  the  leas  t  mrnben  o&  an.cs  between  the  soun.ce  and  the 
sink.  1£  this  chain  has  r  an.cs  then 


max  Q 


_B_ 

rk 


* 


whene  Is  the  capacity  o^  &n  anc  on  Cm,  and  all  othen 
an.cs  have  zeno  capacity. 


11 


3.  FINDING  AN  OPTIMAL  SOLUTION  WHEN  M  0 

While  the  decomposition  approach  could  be  used  to  solve  the 
problem  when  j*  0,  it  would  be  very  inefficient.  Fulkerson's 

modified  primal-dual  algorithm  is  probably  the  most  efficient.  How¬ 
ever  .  an  interesting  alternative  scheme  is  suggested  by  the  approaches 
of  McMasters  and  Mustin  (1970)  and  Doulliez  and  Rao  (1971)  for  looking 
at  problems  involving  capacity  reduction  and  expansion.  These 
approaches  require  the  use  of  a  topological  dual  and,  as  a  conse¬ 
quence,  are  not  particularly  efficient.  They  do,  however,  have  a 
conceptual  appeal. 

Both  papers  present  solution  algorithms  which  involve  dual 
shortest  route  problems  for  identifying  the  cut  sets  of  the  primal 
flow  problem.  They  are,  unfortunately,  restricted  to  planar  net¬ 
works.  This  restriction  has  recently  been  overcome  by  McMasters 
(1971)  who  defines  a  pseudo  topological  dual  associated  with  a  two- 
dimensional  representation  of  the  primal  network.  In  that  paper  the 
special  form  of  the  dual  shortest  route  problem  is  also  stated  and 
an  algorithm  for  solving  the  problem  is  presented.  This  information 
for  both  undirected  and  directed  primal  networks  is  contained  in 
Appendix  A  of  this  paper  for  the  convenience  of  the  reader. 

The  solution  algorithm  presented  below  begins  with  the  con¬ 
struction  of  the  topological  dual  or  psuedo  dual  for  the  primal 

♦ 

network.  The  length  of  the  shortest  and  second  shortest  routes 
through  the  dual  are  then  determined  for  dual  arc  lengths  equal  to 
the  free  capacity  M^  of  the  intersected  primal  arcs.  The  length 


12 


of  the  shortest  route  is  then  increased  by  spending  some  of  the 
resource  budget  on  that  arc  of  the  shortest  route.  If  the  entire 
budget  can  be  spent  without  increasing  the  route's  length  to  that 
of  the  second  shortest  route  then  the  procedure  terminates.  If, 
on  the  other  hand,  there  is  still  a  portion  of  the  budget  left  when 
the  shortest  route  reaches  the  length  of  the  second  shortest  route 
then  the  length  of  the  third  shortest  route  is  determined  and  both 
the  minimum  a^  arcs  of  the  first  and  second  shortest  routes  are 
increased  until  the  budget  is  used  up  or  the  two  routes  attain  the 
length  of  the  third  shortest  route.  The  process  continues  in  this 
manner  until  the  budget  is  used  up. 


Algorithm: 

1.  Assign  the  values  ns  the  arc  flow  capacities  in  the  flow 

network  (primal  network) .  Construct  its  topological  dual  using 
the  procedures  described  in  Appendix  A.  Find  the  shortest  route 
through  the  dual  and  its  length  L^.  Set  n  -  2. 

2.  Find  the  n*"^  shortest  loopless  route  through  the  dual  network 
of  step  1  and  its  length  l/n^ .  Appendix  B  contains  an  n^ 
shortest  loopless  route  algorithm  due  to  Pollack  (1969)  and 
describes  the  modifications  needed  if  the  primal  network  is 
nonplanar . 


n-1 


.(n-1)  _ 


B-  i  H[l(H)-lW] 


k-1 


n-1 

L\ 


k-1 


3.  Compute 
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where  a.  ■  min  a,,  associated  with  the  arcs  on  the  k 

k  ij 


th 


th 


shortest  route  of  the  dual.  In  the  case  where  the  (k-1)  and 
k1"*1  shortest  routes  have  a  "min  a^"  arc  common  then  set 
a^  ■  0.  If  modifications  were  required  in  Pollack's  algorithm 
because  the  primal  network  is  nonplanar  then  any  shortened  arcs 
appearing  on  the  k^  route  should  not  be  included  in  the  a^ 
determination  (see  Appendix  C) . 

4.  Increase  the  lengths  of  the  arcs  associated  Wi  ^,h  . . .  »an_2 

by  the  amount  AY^n^  where 


AY(n)  -  min{L^  -  L(n_1),  B(n_1)}. 


a)  If  AY^  ■  B^n  ^  then  terminate.  The  budget  has  been 
consumed  and  the  optimal  capacity  increase  of  the  primal 
arc  associated  with  af  (r  -  l,2,...,n-l)  is 


*<r)  -  l  « 

J  k-r 


(k) 


For  the  special  case  where  a^  was  set  equal  to  zero  in 
step  3  do  not  compute  Y^  .  All  other  primal  arcs  will 
have  no  increase  in  arc  capacity.  The  value  of  max  Q  is 
equal  to  L(n_1)  +  AY(n). 


b)  If  AY 


(n)  _  L(n)  _  L(n-1)  <  B(n-1)  then  increa8e  n  by 


one  and  return  to  step  2. 


It 


4 .  AN  EXAMPLE 


Consider  the  flow  network  of  Figure  1  with  directed  arcs  and 
nodes  numbered  from  1  to  5.  The  source  is  node  1  and  the  sink 
is  node  5.  The  numbers  on  each  arc  represent  ,  a^j  •  The  v®lue 

of  B  is  assumed  to  be  8.  The  preliminary  steps  of  the  dual  con¬ 
struction  are  shown  on  Figure  1  by  the  dashed  lines  and  nodes  with 
letter  labels. 


SOURCE 


Figure  1. 


The  topological  dual  is  shown  in  Figure  2.  The  numbers  on 
each  arc  represent  the  arc's  length  and  its  associated  a^  value 
in  that  order.  The  dual  origin  is  node  a  and  the  destination  is 
node  e. 

The  shortest  route  through  the  dual  is  a  -  b  -  e  with  a 
total  length  -  3  which  is  the  maximal  feasible  flow  that  the 

primal  network  can  handle  without  paying  for  capacity.  The  second 
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shortest  route  without  loops  is  s-d-c-b-e  with  a  length 
L(2>  - 
(a,b) . 


5.  The  value  of  is  1  and  is  associated  with  arc 


Figure  2. 

Note  that  B(1)  -  -  8.  Because  AY(1)  -  L(2)  -  L(1)  -  2  <  B(1\ 

t:fo  units  of  budget  resource  can  be  consumed  before  the  length  of  the 

(2) 

shortest  route  equals  L  . 

The  third  shortest  route  is  a-d-c-e  with  a  length 

(3) 

L  -  6.  The  value  of  a^i  associated  with  arc  (d,c),  is  2 

(2) 

and  B  ■  2.  We  lengthen  the  arcs  associated  with  a^  and  a ^ 
by  the  amount  AY(2)  -  -  L<2)  -  1. 

(4) 

The  fourth  shortest  route  is  a  -  d  -  e  with  L  ■  8.  We 

get  a^  -  a^  because  arc  (d,c)  has  the  minimum  a^  value  on  both 
(3)  (3) 

routes  and  Bv  '  ■  1.  This  results  in  t Y  -  1  and  the  algorithm 


terminates . 
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The  lengthened  dual  arcs  are  (a,b)  and  (d,c).  Arc  (a,b) 

has  been  lengthened  by  the  amount  AY^  +  AY  ^  +  AY^  -  4.  Arc 

(2)  (3) 

(d,c)  has  been  lengthened  by  the  amount  AY  +  AY  ■  2.  The 
primal  arcs  Intersected  by  (a,b)  and  (d,c)  are  arcs  (1.3)  and 
(3,4)  respectively.  The  capacities  of  these  arcs  have  therefore 
been  Increased  by  4  and  2  respectively.  The  maximal  flow  that 
the  network  can  handle  has  been  Increased  from  3  to  7  units. 

The  optimal  flow  through  each  arc  which  would  give  this  maxi¬ 
mal  flow  can  be  easily  obtained  by  setting  the  arc  lengths  of  (a,b) 
and  (d,c)  at  their  final  values  and  determining  the  shortest  dis¬ 
tances  from  the  origin  to  all  nodes  in  the  dual.  It  has  been  shown 
(Sakarovltch  (1970)),  McMasters  (1971))  that  the  optimal  flow  through 
a  primal  arc  In  a  maximal  flow  problem  will  be  equal  to  the  difference 
between  the  shortest  route  distances  to  the  dual  nodes  incident  with 
the  intersecting  dual  arc.  To  illustrate  this  property  we  consider 
Figure  3  where  arcs  (a,b)  and  (d,c)  have  been  increased  in  length 
according  to  the  optimal  solution  above.  The  value  associated 

with  each  node  of  Figure  3  Is  its  shortest  route  distance  from  node 
a.  Therefore,  the  optimal  flow  through  primal  arc  (.1,3)  is  5 
units  since  -  5.  The  other  optimal  arc  flows  are  shown  in 

Figure  4.  The  numbers  on  each  arc  correspond  to  X^/(M^+Y^j) . 
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5.  EXTENSIONS  OF  THE  PROBLEM 


Nonzero  Lower  Bounds  —  Consideration  of  nonzero  lower  bounds 
on  arc  flow  were  not  explicit  in  Fulkerson's  problem  statement. 
Problems  having  negative  lower  bounds  can,  of  course,  be  handled  by 
replacing  the  primal  arc  by  two  oppositely  directed  arcs  having  zero 
lower  bounds  and  positive  upper  bounds. 

The  general  expression  for  a  lower  bound  constraint  for  the 
parametric  budget  problem  could  be  written  as 


Lu  *  hi +  hi 


(17a) 


where  represents  the  amount  of  reduction  in  L^ 


unit  of  reduction  costs  b 


1J 


If  each 

then  (Id)  would  take  on  the  following 


form: 


1 1*  >Vu  +  Vu1 1 B- 


(17b) 


The  lower  bound  addition  to  the  problem,  regardless  of  the 


sign  of  Ljj  ,  is  easily  handled  by  the  algorithm  since  there  is  always 
a  dual  arc  corresponding  to  the  lower  bound  on  flow  (see  Appendix  A) . 

In  the  example  these  arcs  had  zero  length  since  L^  -  0  for  all  arcs. 

The  b^j  values  are  assigned  to  the  lower  bound  arcs  in  the 
same  way  as  the  a^  values  are  to  the  upper  bound  arcs  in  the 
Appendix  A  procedure.  The  algorithm  then  considers  the  b^'s  as 
merely  a^  values  associated  with  certain  dual  arcs. 

To  prevent  an  arc  associated  with  some  L^  or  from 

being  selected  as  an  a^  arc,  merely  assign  it  a  b^  or  a^ 
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value  very  large.  If  it  is  not  selected  it  will  not  be  changed. 

This  is,  in  fact,  the  reason  for  specifying  infinite  values  of 
a^j  for  the  zero  lower  bound  arcs  in  the  procedures  of  Appendix  A. 

Nonlinear  Convex  a(Y^)  —  When  the  cost  to  increase  capacity 
a(Y^j)  is  nonlinear  but  convex  and  is  zero  for  Y  ■  0  then  a 
piecewise  linear  fit  to  the  cost  function  for  each  arc  can  be  made 
such  that  the  algorithm  can  be  used  with  minor  modifications.  Hu 
(1966)  ,  for  example,  suggests  a  linear  segment  for  each  unit  of 
capacity  change.  Let  a(Y^)  represent  the  general  form  of  the 
cost  function  and  suppose  it  looks  xlke  Figure  5. 


i(r,|) 


Figure  5. 


The  linear  segment  spanning  0  £  Y^  £  y^  would  have  a  slope 


a(y,) 

‘  IT  ' 


the  segment  spanning  Y 2  would  have  a  slope 
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W 


a(y2)  -  a(y1) 
y2  "  yl 


and  so  on. 


The  algorithm  should  begin  by  associating  a^j (y^)  values 
with  each  dual  arc.  It  would  proceed  without  change  until  a  "min 
a^"  arc  reaches  its  y^  value.  If  enough  budget  resource  was 
available  to  further  increase  network  capacity  then  the  "min  a  " 
arc  should  have  its  a^  value  changed  to  a^  ( y and  a  reappraisal 
of  the  a^j  values  for  arcs  of  the  dual  route  currently  under  con¬ 
sideration  should  be  made  to  see  if  some  other  arc  now  has  a  lower 
a^  value.  If  so,  the  length  of  the  new  "min  a^"  arc  should 


be  increased  and  the  first  one  ignored.  The  equation  for  B 


(n-1) 


in  step  3  of  the  algorithm  should  be  modified  to  incorporate  these 
changes  in  a^  and  an  additional  term  should  be  added  to  the 
equation  in  step  4  to  accommodate  the  break  points  y^,  y^t  etc., 
in  the  piecewise  linear  fits  to  the  a(Y^)  curves  of  the  dual  arcs 
being  lengthened. 

These  modifications  would  be  directly  applicable  to  problems 
having  linear  a(Y^)  cuLves  but  upper  bounds  on  Y^  values. 


An  Interdiction  Problem  —  The  interdiction  problem  studied 
by  McMasters  and  Mustin  (1970)  seeks  to  spend  money  to  reduce  the 
maximum  possible  flow  of  an  enemy's  supply  network,  ’his  is  done 
by  reducing  the  capacities  of  certain  arcs  of  the  network.  If  Z^j 
is  the  amount  of  capacity  reduction  of  arc  (i,j)  and  a^  is  the 
unit  reduction  cost  then  (Id)  would  have  the  form 
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*  k  V«  *  B- 


If  Is  the  arc  capacity  before  interdiction  and  i  0  is 

the  least  possible  capacity  after  interdiction  then  (18)  and  (19) 


describe  the  bounds  on  and  . 

0  S  Zij  *  Mij  "  ®ij  * 

(19) 

0  *  Xij  *  Mij  ‘  Zij  • 

(20) 

The  algorithm  of  section  3  is  easily  modified  to  solve  this 
problem.  The  a^  arc  of  the  dual  shortest  route  would  be  shortened 
in  length  until  it  reached  its  m^  value  or  the  budget  is  spent. 

If  the  m^  value  is  reached  first  then  shorten  the  arc  having  the 
next  smallest  value  until  it  attained  its  m^  value  or  the 

budget  is  spent.  Repeat  this  process  for  each  successive  arc  of  the 
dual  shortest  route.  The  final  length  of  the  f^ute  should  then  be 
recorded . 

The  dual  second  shortest  route  should  be  analyzed  in  the  same 
fashion  and  its  finsl  length  compared  with  that  of  the  shortest  route. 
That  route  having  the  shortest  final  length  is  retained  for  further 
comparisons.  This  process  must  be  repeated  for  each  route  through 
the  dual.  The  shortest  final  route  crosses  the  primal  arcs  to  be 
interdicted.  The  amount  of  effort  allocated  to  each  primal  arc  depends 
on  the  length  of  the  intersecting  dual  arc. 
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The  algorithm  presented  in  the  McMaster3  and  Mustin  paper  has 
this  same  flavor  but  is  more  efficient  if  ra^  >  0  for  some  arcs. 

Thel*-  procedure  begins  by  finding  the  shortest  route  through  the 
dual  when  all  arcs  of  the  dual  are  set  at  their  m^  values.  The 
amount  of  money  needed  to  attain  the  shortest  route  length  is  then 
computed.  If  it  exceeds  the  budget  then  some  "unspending"  is  required. 
If  it  does  not  exceed  the  budget  the  problem  is  solved. 

The  "unspending"  looks  for  the  most  expensive  arc  (that  corre¬ 
sponding  to  max  a^j)  on  the  dual  shortest  route.  The  length  of 
that  arc  is  then  increased  until  it  reaches  its  M  value  or  the 

budget  constraint  is  satisfied.  If  it  reaches  its  valuci  first 

then  the  next  most  expensive  arc  is  also  lengthened.  The  process 
repeats  until  the  budget  constraint  is  satisfied.  The  firal  length 
of  the  shortest  route  is  recorded. 

The  second  shortest  dual  route  based  on  m^  values  is  next 
determined  and  its  length  checked  against  the  final  length  of  the 
shortest  route.  If  it  is  longer  than  that  final  length  the  problem 
is  solved;  otherwise  one  or  more  arcs  of  the  second  shortest  route 
are  lengthened  to  meet  the  budget.  The  final  length  of  the  second 
shortest  route  is  compared  with  that  of  the  shortest  and  the  route 
with  minimum  length  is  retained  for  further  comparisons.  The  process 
continues  until  all  dual  routes  have  been  examined  or  some  dual  route 
having  all  arcri  at  their  m^  values  exceeds  in  length  the  shortest 
preceding  route  meeting  the  budget  constraint. 
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6 .  APPENDIX  A 

Construction  of  the  Topological  Dual  —  The  original  flow  network 
will  be  called  the  primal  network.  A  mesh  of  a  planar  primal  net¬ 
work  is  any  region  surrounded  by  nodes  and  arcs  but  containing 
neith''-:  ip  the  plane  on  which  the  network  is  constructed.  The  region 
of  the  plane  completely  jjTounuing  the  primal  network  will  be  called 
the  external  me-'h.  The  construction  of  the  dual  of  a  source-sink 
planar  directed  network  consists  of  the  following  steps  (McMasters 
(1971)). 

1.  Denote  the  original  maximal  flow  network  as  the  primal  network. 
Connect  an  artificial  arc  between  the  sink  and  source  of  the 
primal  and  position  it  below  the  network.  The  resulting  network 
will  be  referred  to  as  the  modified  primal  network. 

2.  Place  a  node  in  each  mesh  of  the  modified  primal  including  the 
external  mesh.  Let  the  origin  of  the  dual  be  the  node  in  the 
mesh  involving  the  artificial  arc  and  the  destination  be  the  node 
in  the  external  mesh. 

3.  For  each  arc  in  the  primal  (except  the  artificial  arc)  construct 
two  oppositely  directed  arcs  that  intersects  it  and  joir.  with 
nodes  in  the  meshes  adjacent  to  it. 

4.  Assign  the  value  of  the  upper  bound  capacity  of  the  primal  arc 
as  the  length  of  the  intersecting  dual  arc  having  the  same 
direction  that  the  primal  arc  would  have  if  it  were  rotated  90° 
counterclockwise.  Assign  to  the  oppositely  directed  dual  arc  a 
length  equal  to  the  negative  of  the  lower  bound  capacity  of  the 
primal  arc  (these  lengths  will  all  be  zero  for  problem  (1)). 
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If  the  primal  network  i9  not  source-sink  planar  then  modify 
step  1  above  to  read: 

1.  Construct  a  two-dimensional  representation  of  the  flow  network 
such  that  all  arcs  are  straight  lines.  Assign  a  psuedo  node  to 
every  intersection  of  arcs  not  at  a  node  in  this  representation. 
Connect  an  artificial  arc  between  the  sink  and  source  of  this 
representation  and  position  it  below  the  network.  The  resulting 
network  will  be  referred  to  as  the  modified  prljnal  network. 

The  remainder  of  the  steps  for  dual  construction  are  the  same  as 
above. 


Assignment  of  a^  Value  — 
primal  network  to  those  dual 
values.  Assign  a^  ■  «  to 


Assign  the  a^  values  of  the  original 
arcs  having  lengths  equal  to  the  M 
all  dual  arcs  corresponding  to  the  lower 


bound  capacity  of  the  primal  arc. 
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7.  APPENDIX  B 

An  nC^  Shortest  Looplesa  Route  Algorithm  —  We  are  interested  in 
looplcss  routes  (contain  no  cycles)  through  the  dual  because  we  want 
those  routes  corresponding  to  primal  cut  sets^  which  disconnect  the 
primal  into  two  subgraphs,  one  containing  the  source  *nd  the  other 
the  sink.  We  know  that  any  route  containing  a  cycle  would  not  corre¬ 
spond  to  such  a  cut  set  (Ford  and  Fulkerson  (1956)).  We  are  also 
interested  in  having  an  algorithm  which  allows  us  to  get  the  next 
shortest  route  at  any  time.  Pollack  (1969)  developed  the  following 
algorithm  which  has  both  of  these  feature;.. 

1.  Determine  the  shortest  rcute  through  the  network  using  an  algo¬ 
rithm  such  as  Dijkstra's  if  all  arc  lengths  are  non-negative  or 

2 

Yen's  if  some  arc  lengths  are  negative  (see  Dreyfus  (1969)).“ 

2.  To  determine  the  second  shortest  route  remove  the  first  arc  from 
the  shortest  route  and  solve  for  the  shortest  route  through  the 
remaining  network.  This  route  is  a  candidate  for  the  second 
shortest  route.  Record  the  route  and  its  length  and  replace  the 
arc.  Remove  the  second  arc  of  the  shortest  route,  solve  for  the 
shortest  route  through  the  remaining  network,  record  its  length 
and  replace  the  second  arc.  Continue  this  process  until  all  arcs 

^  A  cut  set  of  a  connected  graph  is  defined  to  be  a  disconnecting  set 
of  arcs  which  contains  no  proper  subset  which  also  disconnects  the 
graph. 

2 

Negative  arc  lengths  will  occur  in  the  dual  networks  of  problems 
having  positive  lower  bounds  on  arc  flows.  An  infeasible  flow  is 
detected  by  a  cycle  of  negative  length  in  the  dual  (McMasters  (1971)). 
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on  the  shortest  route  have  been  removed  and  replaced.  Examine 
the  list  of  second  shortest  route  candidates.  That  candidate 
having  the  shortest  length  Is  the  second  shortest  route.  In  the 
case  of  a  tie  between  two  candidates  arbitrarily  select  one  as 
the  second  shortest  and  specify  the  other  as  the  third  shortest. 

3.  If  there  were  no  ties  for  the  second  shortest  route,  the  third 
shortest  route  candidates  are  determined  by  removing  an  arc  from 
the  shortest  and  an  arc  from  the  second  shortest  routes  and  solving 
for  the  shoicesi.  route  through  the  remaining  network.  This  must 

be  done  for  all  combinations  of  two  arcs  from  the  shortest  and 

second  shortest  routes.  The  list  of  candidates  is  then  examined 

and  the  third  shortest  route  Is  that  candidate  having  the  shortest 

length. 

th 

4.  The  n  shortest  route  is  obtained  by  first  finding  the  shortest 
route  through  the  remaining  network  after  each  combination  of 
arcs,  one  from  each  of  the  preceding  n-1  shortest  routes,  has 
been  removed.  That  candidate  having  the  shortest  length  is  the 
nfck  shortest  route. 

Modifications  for  a  Nonplanar  Primal  —  If  the  primal  netwerk  is  not 
source-sink  planar  then  loopless  routes  through  the  dual  may  not 
correspond  to  cut  sets  of  the  original  primal.  Modifications  to 
Pollack's  algorithm  may  therefore  be  required  to  insure  selection 


of  the  correct  dual  routes. 
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Let  the  dual  nodes  be  numbered  1,2,...,N  where  1 
corresponds  to  the  dual  origin  and  N  corresponds  to  the  dual 
destination  or  terminal  node.  Let  the  group  of  four  dual  nodes 
surrounding  a  primal  pseudo  node  ^  be  represented  by  (i,J,k,m)^. 
Let  ?  be  the  set  of  all  such  groups  in  a  dual  network.  Let 
represent  the  permanent  label  on  dual  node  i  after  the  shortest 
route  from  the  dual  origin  (V^  -  0)  to  every  dual  node  has  been 
determined.  Let  be  the  length  of  the  dual  arc  (i,J).  If, 

during  each  step  of  Pollack's  algorithm, 

vj  -  vi  ■  v=  -  \  <18) 

for  all  node  sets  of  9  then  no  modifications  are  necessary.  If 
(18)  is  not  satisfied  for  some  set  then  use  the  following  subroutine 
to  change  the  values.  McMasters  (1971)  explains  the  reasoning 

behind  condition  (18)  and  proves  that  this  subroutine  does  accomplish 
the  desired  result. 

Adjustment  Subroutine: 

1.  Compute  a  set  of  numbers  Y.  (i  •  1,2,...,N).  where  Y„  •  Vu 

1  N  N 

and 

\  •  “*<W 

for  i  ■  1,2, . . . ,N-1. 
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2.  For  each  set  (i,J,k,ra)  for  which  (18)  Is  not  satisfied  compute 

9 

the  following  numbers: 


■  "ln(,ij'vJ-',i'VYk) 

l2U)  - 


If  1^(9)  <  i2(9)  ^or  one  or  more  sets  then  go  to  step  3; 
otherwise,  go  to  step  A. 

3.  Determine  that  set  (i,j,k,m)*  corresponding  to 

♦ 


12U)  -  ^(9) 

and  charge  , 
with  that  set  to 


-  maxU2(4i)  -  i1(v)  I  ^(9)  <  ^2(^)) 


and  l  .  of  tl—  arcs  associated 
mk 


^(9) 


ji 


Go  to  step  5. 

A.  Select  any  set  (i,j,k,m)  ,  for  which  (18)  is  not  satisfied  and 

9 

change  the  values  of  l..,  1..,  I,  ,  and  l  .  on  the  arcs 

ij  ji  km  mk 

associated  with  that  set  to 


^(9), 


ji 


-^(9). 
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5.  Recompute  the  shortest  routes  from  the  dual  origin  to  each  node. 

(a)  If  (18)  is  not  satisfied,  return  to  step  1  of  this  subroutine. 

(b)  If  (18)  is  satisfied,  record  the  shortest  route  through  the 

dual  and  its  length.  Identify  any  shortened  arcs  on  this 
route.  If  this  route  is  later  found  to  be  the  shortest 

route  then  these  arcs  should  not  be  "removed"  on  subsequent 
iterations  of  Pollack's  alforitha.  Return  all  shortened 
arcs  to  their  original  lengths  and  continue  with  Pollack's 
algorithm. 

The  adjustment  subroutine  artificially  reduces  the  lengths 
of  certain  arcs  of  the  dual.  At  the  end  sf  the  subroutine  one  or 
more  of  these  arcs  may  appear  on  the  shortest  route.  In  searching 
for  the  next  shortest  route  using  Pollack's  algorithm,  these 
shortened  arcs  should  not  be  removed  since  "removal"  is  normally 
accomplished  in  shortest  route  algorithms  by  assigning  an  infinite 
length  to  the  arc.  Increasing  the  original  length  of  a  shortened 
arc  would  result  in  the  route  just  determined  by  the  adjustment  sub¬ 
routine  being  also  determined  as  a  candidate  for  the  next  shortest 


route. 
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8.  APPENDIX  C 

Selection  of  a^  Area  for  a  Nonplanar  Primal  —  The  adjustment 
subroutine  of  Appendix  B  may  cause  arcs  on  a  dual  kC^  shortest 
route  to  be  shortened.  Because  increasing  the  original  length  of 
a  shortened  arc  would  not  increase  the  length  of  the  kC^  shortest 
route  we  would  gain  no  increased  primal  flow  capacity.  Therefore, 
we  should  select  that  arc  corresponding  to  "min  a^"  from  only 
those  arcs  on  the  route  which  were  not  shortened  by  the  subroutine. 
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